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Q-ball is a non-topological soliton whose stability is ensured by global U(l) symmetry. We 
study a Q-ball which arises in the Affleck-Dine mechanism for baryogenesis and consider its possible 
instability due to 1/(1) breaking term (^4-term) indispensable for successful baryogenesis. It is found 
that the instability destroys the Q-ball if its growth rate exceeds inverse of the typical relaxation 
time scale of the Q-ball. However, the instability is not so strong as it obstructs the cosmological 
formation of the Q-balls. 



I. INTRODUCTION 

The Q-ball is a non-topological soliton that arises in 
scalar field theory with some global U(l) symmetry 0. 
The Q-ball solution naturally exists in the minimal super- 
symmetric (SUSY) standard model (MSSM) 0HI1E1, 
especially in the context of the Affleck-Dine (AD) baryo- 
genesis |6j where the MSSM flat directions play impor- 
tant roles in baryon number generation. In this case the 
Q-ball consists of squarks and sleptons, therefore carry- 
ing baryonic and/or leptonic charges. It was shown that 
Q-balls with large baryon number are actually produced 
in the early universe in Refs. @, @- Furthermore, a Q- 
ball with large baryon number is stable against decay into 
protons in theories based on the gauge-mediated SUSY 
breaking. Therefore, such stable Q-ball can be a promis- 
ing candidate for dark matter as well as the source of the 
baryon number of the universe, which makes the Q-balls 
very attractive 0, 0, 0|. Q-balls in the gravity- 
mediated SUSY breaking models, on the other hand, are 
unstable and decay into the standard particles, producing 
the lightest SUSY particles (LSPs). Then it is possible 
that the Q-balls account for the dark matter (=LSPs) 
and the baryon asymmetry of the universe simultane- 
ously la . For further applications of Q-balls and their 
variant [13, see e.g. Refs. [II El El Ej . 

In the AD mechanism we need the U(\)b(l) breaking 
terms for successful generation of the baryon (lepton) 
number. Here U(1)b(l) is the global symmetry associ- 
ated with baryon (lepton) number, and in the follow- 
ing we drop the subscript B(L), since the distinction 
makes no difference to the following discussion. If Q- 
balls are not formed, the U(l) breaking terms can be ne- 
glected soon after the baryon number is generated. This 
is because the cosmic expansion decreases the amplitude 
of the AD field and the U(l) breaking terms become 
much smaller than the U(l) conserving ones. However, 
once the Q-balls are formed, since the amplitude of the 
AD field inside the Q-ball is fixed it becomes nontrivial 
whether the U(l) breaking terms are truly negligible or 



not. Therefore, in the present paper, we study the pos- 
sible instability due to the £7(1) breaking A-terms and 
their effect on the Q-ball stability. It is found that the 
A-term induces instability similar to that of the paramet- 
ric resonance and it upsets the stability of the Q-ball if 
the growth rate of the induced instability is larger than 
the inverse of the relaxation time of the Q-ball configu- 
ration. However, in the realistic cosmological situation, 
the A-term inside the Q-ball is not so strong to cause 
the strong instability, hence the previous studies which 
assumed the Q-ball stability remain valid. 

II. LINEAR ANALYSIS ON INSTABILITIES 

First let us consider the instabilities due to U(l) break- 
ing term in the homogeneous background. To this end 
we perform linear analysis assuming small perturbations. 
The potential of the AD field <& is written as 

Vm = m|(l + inog(ff))|4>| 2 

where m$ is the mass of the AD field, K a numerical 
coefficient of the one- loop correction, m 3 / 2 the gravitino 
mass, M the renormalization scale, and M* some cut-off 
scale for the nonrenormalizable operator. Also A is as- 
sumed to be a O(l) real parameter, and K is assumed to 
be negative. One can see that the second term (called 
A-term) breaks the U(l) symmetry. In the following 
argument we can safely neglect the last term since the 
amplitude of the AD filed is relatively small. 

In order to study the instabilities of the AD field we 
first divide the AD field into homogeneous part and fluc- 
tuation: <f> = $ + 5$. The equations of motion are 

# + 3H$ + m% (l + K + if log f^a )) * 

+ ^3/2-^ = 0, (2) 



for the homogeneous mode, and 

S$ + 3H8& + -^5$ + Ami$ 



+ m% 1 + K + Klog 
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(3) 



for the fluctuation with the wavenumber k in the momen- 
tum space. Although we have introduced the Hubble pa- 
rameter H and the scale factor a in the above equations, 
we will neglect the cosmic expansion for the moment. 

In order to parametrize the strength of the instability, 
let us define £ as the ratio of the A-term to the mass 
term: 



d-2 



dm%M£- 3 



(4) 



where $o is the maximal magnitude of $ during the 
course of the oscillation in the limit of £ — > and we 
will set M = $o. Note that the amplitude of the oscil- 
lation in the limit of £ — * is constant if we neglect the 
cosmic expansion. We have numerically solved Eqs. (J5J) 
and © by decomposing $ and <5$ into their real and 
imaginary components: 



<5$ 



)(f>i + iS<j>2 
V2 ' 



(5) 



For convenience we also define the polar decomposition: 
$ = <j)e 10 /y/2. We have normalized the dimensionless 
variables as ifi = 0i/m$, 5<fi = 5<f)i/m$, (i = 1,2), k = 
h/m®, t = m$>t, and Xi = m<s>Xi (i = 1> • • • , D), where 
D is the spatial dimension. Furthermore we set K = 
in order to concentrate on the instability that originates 
from the A-term. Fig^ shows the instability bands for 
d = 4 and £ = 0.1 with the initial condition 9 = ir/8 
and 9' = 1, where the prime denotes the differentiation 
with respect to r. Notice that with the above initial 
values the orbit of the homogeneous mode is close to 
a circle in the (cfrijfa) plane. From the figure one can 
see that there exist instability bands similar to those of 
the parametric resonance; the strongest instability band 
is located at k — ~ 0.8 and several weaker ones at 
shorter wavelengths. We also show the result for d = 6 
and £ — 0.1 with the initial condition 6 = ir/12 and 
9' = 1 in Fig[2 In this case the strongest instability is 
at larger k around 1.7. 

In Figs. |2| and 0] the instability bands are shown for 
the elliptic orbit of the AD field. We have chosen the 
initial conditions as 6 = ir/8 and 9' = 0.01 for d = 4 




FIG. 1: Instability bands due to the A-term for d = 4, £ = 0.1 
and circular orbit. 




FIG. 2: Same as Fig. □ but d = 6. 



and 9 = 7r/12 and 9' = 0.01 for d — 6, respectively. One 
can see that the instability bands are located at different 
wavenumber, and that the instabilities are weaker and 
bands narrower, compared to the cases of the circular 
orbit. 

In order to express the strength of the strongest in- 
stability in the case of the circular orbit, we introduce 
the growth rate T g ; the fluctuation at the peak of the 
strongest instability band grows as oc exp(r g i). The 
growth rate depends on £ and we find that r g is pro- 
portional to £ as shown in FigO Thus we obtain 



(6) 



with 74 w 2.0 and « 3.5. 

Before closing this section, let us consider the effect 
of the cosmic expansion on the instabilities. When the 
oscillation of the AD field starts, the Hubble parameter 
H is comparable to m$. Hence the growth rate of the 
A-term instability is comparable to H only if £ ~ O(l). 
However, since the amplitude of the AD field decreases as 
a -3 / 2 , £ quickly becomes smaller (oc a~ 3 ( d ~ 2 )/ 2 ) and the 
growth rate is soon overcome by the cosmic expansion. 
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FIG. 3: Same as Fig.Qbut the elliptic orbit. 




FIG. 4: Same as Fig.EJbut the elliptic orbit. 



Therefore, it seems that the A-term instability does not 
play an important role at least in the Q-ball formation. 
However, it is still nontrivial how the instability affects 
the evolution of Q-balls after they are formed, since the 
amplitude inside them is fixed. 



III. INSTABILITIES INSIDE Q-BALL 

In the previous section we have found that the A-term 
causes the instabilities in the homogeneous motion of the 
AD field. Next let us investigate the instabilities due to 
the A-term inside the Q-ball. Without the A-term, the 
Q-ball solution for the potential JIJ is obtained with use 
of the Gaussian ansatz : 



0(r) = <£(0)e~^r, 



(7) 

(8) 
(9) 



FIG. 5: ^-dependence of the growth rate of the most amplified 
mode, F 9 . The +'s and x 's are the numerical results for d — 4 
and d = 6, respectively. The solid line denotes r 9 /m$ ~ 3.5£, 
while the dotted one denotes r 9 /m$ ~ 2.0£. 



where the Q-ball radius R and angular velocity u> are 
given by 



R- 



m%\K\ 



m\{l + 2\K\ 



(10) 
(11) 



To study the instabilities due to the A-term inside the Q- 
ball, we have numerically solved the evolution of the AD 
field on one and two dimensional lattices. The equation 
of motion for the potential is 



$ - V 2 $ + m 2 



Am 



3/2 mF* 



1 



o. 



(12) 



We have neglected the expansion of the universe since the 
Q-balls decouple from the cosmic expansion once they are 
formed. 

In the numerical calculations we have set M = 
</>(0)/V2 and investigated several values of K between 
—0.01 and —0.3 for both d = 4 and d = 6, varying £. For 
the Q-ball configuration, £ is evaluated at the center of 
the Q-ball; <I>o in Eq. (0} should be interpreted as <& c , the 
field value at the center. Also we have added small initial 
seed for fluctuations of O(10~ 5 ). The size of the lattices 
is 2048 and 200 x 200 for ID and 2D, respectively. A 
number of simulations with varying values of the lattice 
spacing and box size have been done to ensure that the 
exact values of the lattice parameters do not affect the 
results. 

Then we have found that the Q-ball breaks up for £ 
larger than O(10" 2 ) (see Tables |T| and HJ for the precise 
values). Figs. El [7| and |S1 show how the Q-balls are dis- 
persed by the A-term instabilities. It is worth noting that 
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FIG. 6: The typical breakdown of the Q-ball due to the A- 
term for d = 4 on 1 + 1 lattices, r — 0, 10 3 and 10 4 from 
top to bottom. The adopted parameters are K = —0.01 and 
i = 8 x 10" 3 . 
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FIG. 7: Same as Fig. but d = 6. r = 0, 3.0 x 10 3 and 
5.0 x 10 from top to bottom. The adopted parameters are 
K = -0.04 and £ = 4.5 x 10~ 2 . 



TABLE I: The critical values £ c for d = 4 



Criterion 


K = -0.01 


K = -0.03 


K = -0.1 


K = -0.3 


A 


4 x 10" 3 


9 x 10" 3 


2.8 x 10" 2 


6.4 x 10" 2 


B 


3 x 10" 3 


9 x 10" 3 


2.6 x 10" 2 


5.5 x 10 -2 


TABLE II: The critical values £ c for d = 6 


Criterion 


K = -0.01 


K = -0.03 


K = -0.1 


K = -0.3 


A 


2.2 x 10~ 2 


3.5 x 10~ 2 


2.3 x 10~ 2 


1.2 x 10~ 2 


B 


2.1 x 10~ 2 


3.5 x 10~ 2 


2.3 x 10~ 2 


1.2 x 10~ 2 



the Q-ball is actually stable for small enough values of £. 
The critical values (£ c ) above which the Q-ball becomes 
unstable due to the A-term instabilities for several val- 
ues of K are shown in Tables |U and |n] We have followed 
the evolution until r = 10 4 and then decided whether 
the Q-ball configuration is lost or not on the basis of the 
following criterion; if the total charge of the Q-ball is less 
than 10% (50%) of the initial value, we judge that the 
Q-ball is dispersed. We call this criterion A (B). 

The critical value £ c is found to depend on the values 



of both d and K to & certain degree. The dependence of 
£ c on d comes from the fact the instability bands have 
different structures between the cases of d = 4 and d = 6 
(see Figs, n and El- The dependence on the other vari- 
able, K, possibly comes from the following two effects; 
(i) the Q-ball is more stable for larger \K\; (ii) the Q-ball 
is more tolerant to perturbations if its charge is larger, 
and the charge becomes larger for smaller \K\ with the 
amplitude 4>(0) fixed because of the dependence of the ra- 
dius on \K\ (see Eq. IjlOfl L Such a naive reasoning might 
explain the behavior of £ c especially in Table UTI 

Apart from the detailed dependence, can we tell a 
rough value of £ c ? Put another way, is £ c ~ O(10 -2 ) 
reasonable? It is conceivable that £ c is determined by 
the competition between the growth of the instabilities 
and the relaxation of the Q-ball configuration. Since the 
Q-ball configuration minimizes the energy of the system 
in the limit of £ — > 0, the Q-ball tends to keep its config- 
uration for a certain range of £. It is expected that the 
relaxation time scale of the Q-ball is set by the mass of 
the AD field, m$, possibly multiplied by some powers of 
\K\. On the other hand, the growth rate calculated from 
Eq. 10 is about 0.1m$ for the critical values of £. So, 
the relaxation time scale of the Q-ball should be around 
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FIG. 8: The typical breakdown of the Q-ball due to the A- 
term for d = 6. r = 0, 1.0 x 10 3 and 1.5 x 10 3 from top 
to bottom. The parameters we adopted are K — —0.04 and 
£ = 4.0 x 10" 2 . 



10 . This relaxation time scale is in an agreement 
with Ref. where the excited Q-balls were studied by 
numerical simulations. Therefore it is probable that the 
instability due to the A-term destroy the Q-ball if its 
growth rate exceeds the inverse of the typical relaxation 
time scale of the Q-ball. 

Lastly let us comment on the realistic value of £ at the 
Q-ball formation. When the AD field starts oscillating, 
£ should be order unity if m 3 / 2 — m$ and A ~ 0(1). 
Then the typical value of £ is given by £ ~ {& c /®osc) d ~ 2 , 



where <I? osc is the amplitude of the AD field at the onset of 
oscillation. According to the numerical calculations (see 
e.g. Ref. H3), $ c /$ OS c = 6x 10~ 3 (|if|/0.1) 3 / 4 . Thus 
the realistic value of £ is much smaller than £ c : 



3 x 10" 



10" 



\K\ 
0.1 



\K\ 
0.1 



for d = 4 



for d = 6 



(13) 



Note that £ becomes smaller while d increases. 

IV. CONCLUSION 

In this paper we have investigated how the A-term af- 
fects the evolution of the AD field, especially paying at- 
tention on the stability of the Q-balls. In the linear anal- 
ysis we first have found that there exist instability bands 
similar to those of parametric resonance. The growth 
rate of the instabilities, however, is not so large after the 
baryon asymmetry is generated. Thus the Q-ball forma- 
tion in the expanding universe would not be disturbed by 
the presence of the instabilities. Next we have studied the 
stability of Q-balls and estimated the critical value of £ c , 
above which the Q-balls cannot stay stable, for several 
values of K, From our result it is conceivable that £ c is 
determined by the competition between the growth rate 
and the relaxation rate of the Q-balls. It should be also 
noted that the obtained £ c is rather large in the realistic 
cosmological situations. Therefore the extensive studies 
on Q-balls thus far should remain valid, since the realistic 
value of £ should be smaller than £ c . On the other hand, 
the instability found in the present paper may be im- 
portant when Q-balls grow up by absorbing t/(l) charge 
such as solitosynthesis |l9j . 
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